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An  algorithm  is  described  and  evaluated  for  determining  the  absorption  and  backscattering  coefficients 
a(z)  and  bb(z)  from  measurements  of  the  nadir-viewing  radiance  Lu(z)  and  downward  irradiance  Ed(z). 
The  method,  derived  from  radiative  transfer  theory,  is  similar  to  a  previously  proposed  one  for  Eu(z)  and 
Ed(z),  and  both  methods  are  demonstrated  with  numerical  simulations  and  field  data.  Numerical 
simulations  and  a  sensitivity  analysis  show  that  good  estimates  of  a(z )  and  bb(z)  can  be  obtained  if  the 
assumed  scattering  phase  function  is  approximately  correct.  In  an  experiment  in  Long  Island  Sound, 
estimates  of  a(z)  derived  with  these  methods  agreed  well  with  those  obtained  from  an  in  situ  reflecting 
tube  instrument.  ©  1999  Optical  Society  of  America 
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1.  Introduction 

The  estimation  of  the  beam  absorption  and  backscat¬ 
tering  coefficients  a  and  bb  of  natural  waters  is  a 
primary  goal  of  optical  oceanographers.  These  in¬ 
herent  optical  properties  (IOP’s)  determine  the  ocean 
surface  color,  the  quality  of  underwater  visibility,  and 
the  transfer  rates  of  heat  and  photosynthetically 
available  radiation  into  the  upper  ocean.  The  values 
of  the  IOP’s  depend  on  the  type  and  concentration  of 
the  water  components,  primarily  of  phytoplankton, 
suspended  sediment,  and  dissolved  and  particulate 
organic  material.  Therefore  multispectral  measure¬ 
ments  of  a  and  bb  can  be  used  to  determine  the  water 
constituent  types  and  abundances.  The  values  of  a 
and  bb  can  be  estimated  from  remote-sensing  mea¬ 
surements  of  ocean  color;  however,  in  situ  determi¬ 
nations  are  more  accurate  and  are  required  for  the 
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development  and  verification  of  remote-sensing  algo¬ 
rithms. 

There  are  two  general  approaches  to  in  situ  IOP 
determination.  The  first  is  to  measure  the  attenua¬ 
tion  of  an  active  light  source  through  a  water  sample 
with  an  instrument  that  physically  isolates  the  water 
sample  from  the  surrounding  environment.1-2  The 
most  widely  used  of  these  instruments  is  the  ac-9 
(WetLabs,  Inc.,  Philomath,  Oregon),  which  is  capable 
of  measuring  fine-scale  vertical  variations  in  the 
IOP’s  and  can  operate  day  or  night,  within  or  below 
the  euphotic  zone.  It  has  the  disadvantages, 
though,  that  corrections  must  be  made  to  the  esti¬ 
mates  of  a  to  reduce  scattering  errors  and  that  the 
volume  of  the  sample  is  smaller  than  what  is  desired 
for  most  oceanographic  applications.  Also,  the  small 
intakes  of  these  instruments  can  break  up  or  fail  to 
collect  large  suspended  aggregates  that  affect  under¬ 
water  visibility. 

The  second  approach  to  IOP  determination  is  to 
measure  the  natural  light  within  the  water  column 
and  invert  the  radiative  transfer  equation  to  calcu¬ 
late  the  values  of  a  and  bb  that  are  consistent  with  the 
observed  light  measurements.  This  requires  that 
the  effects  on  the  light  field  of  the  desired  IOP’s  be 
separated  analytically  or  empirically  (rather  than 
physically)  from  the  effects  of  the  illumination,  sea 
state,  and  bottom  conditions.  This  approach  en¬ 
ables  one  to  obtain  the  water  and  light  field  proper¬ 
ties  from  the  same  instrument,  but  it  can  be 
employed  only  during  the  day  and  in  the  upper  layers 
of  the  ocean  where  there  is  sufficient  natural  illumi¬ 
nation.  With  the  algorithm  presented  here,  each 
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IOP  estimate  is  affected  by  an  optically  large  volume 
of  water  surrounding  the  nominal  location  of  the  es¬ 
timate.  Therefore  it  is  not  possible  to  resolve  the 
line  vertical  structure  of  waters  that  have  high  spa¬ 
tial  variability.  Nonetheless,  even  in  highly  vari¬ 
able  waters,  the  method  can  retrieve  spatially 
averaged  estimates  of  the  IOP’s  that  may  be  appro¬ 
priate  for  remote  sensing  and  other  large-scale  appli¬ 
cations. 

Because  the  two  general  approaches  for  IOP  deter¬ 
mination  are  independent  of  one  another,  it  should  be 
possible  to  demonstrate  optical  closure  between  the 
two.  This  is  attempted  here  with  data  collected  in 
Long  Island  Sound  by  comparing  the  estimates  of 
a(z )  obtained  from  an  ac-9  instrument  with  those 
obtained  from  measurements  of  the  downward  irra- 
diance  EJz)  and  either  the  nadir-viewing  radiance 
LJz)  or  the  upward  irradiance  EJz)  at  several 
depths  z. 

The  method  for  determining  a(z )  and  bjz)  from 
profiles  of  EJz)  and  EJz)  has  been  reported  previ¬ 
ously.3  The  approach  consists  of  estimating  the  lo¬ 
cal  asymptotic  values  of  the  downward  diffuse 
attenuation  coefficient  and  the  irradiance  reflectance 
KJz)  and  RJz),  respectively,  and  then  calculating 
a(z)  and  bb{z)  from  KJz)  and  RJz).  The  method 
was  evaluated  with  numerical  simulations  and  a  sen¬ 
sitivity  analysis.  Here  the  algorithm  is  extended  to 
allow  use  of  LJz)  instead  of  EJz)  and  is  applied  to 
both  simulations  and  field  data.  Also  it  is  shown 
that  the  iteration  step  of  the  algorithm  (for  determin¬ 
ing  the  single-scattering  albedo  w0)  can  be  replaced 
by  a  simple  inversion  of  the  one-to-one  correspon¬ 
dence  between  co0  and  R,_  that  exists  for  any  given 
scattering  phase  function. 

Gordon  and  Boynton4-5  have  presented  an  alterna¬ 
tive  method  for  the  estimation  of  a(z)  and  bjz)  from 
Ed(z )  and  EJz)  or  LJz).  With  this  method,  esti¬ 
mates  of  a(z)  and  bjz)  are  used  to  compute  profiles  of 
EJz)  and  Eu(z)  or  LJz)  with  a  radiative  transfer 
numerical  code.  The  differences  between  the  com¬ 
puted  and  the  measured  Ed(z)  and  EJz)  or  Lu{z )  are 
used  to  iteratively  adjust  a(z)  and  bjz)  until  an  op¬ 
timal  solution  is  obtained.  This  approach  is  far 
more  computationally  intensive  than  the  method  de¬ 
scribed  here  because,  with  our  method,  no  iteration  is 
required.  Also  the  method  of  Gordon  and  Boynton 
requires  the  surface  conditions  (illumination  and  sea 
state)  as  an  input,  whereas  our  method  does  not  re¬ 
quire  surface  information.  Although  this  makes  the 
Gordon  and  Boynton  algorithm  more  complicated  to 
apply,  the  additional  information  should  make  it 
more  accurate  than  our  algorithm  at  shallow  optical 
depths.  Both  approaches  can  accommodate  the  ef¬ 
fects  of  bottom  reflectance  in  optically  shallow  wa¬ 
ters. 

Our  Lu-Ed  method  is  presented  in  Sections  2  and  3, 
and  a  sensitivity  analysis  and  numerical  testing  of 
the  method  is  given  in  Sections  4  and  5,  respectively. 
The  application  of  both  the  Eu-Ed  and  Lu-Ed  algo¬ 
rithms  to  field  data  is  demonstrated  in  Section  6,  and 
a  discussion  is  given  in  Section  7. 


2.  Basic  Equations 

We  are  interested  in  the  azimuthally  averaged  radi¬ 
ance  L(z,  |x)  (in  units  of  W  m  2  sr  1  nm  )  that 
satisfies  the  integrodifferential  radiative  transfer 
equation 

p,8L(z,  I ±)/dz  +  c(z)L{z,  |x) 

=  6(z)J  (3(z,  p-,  p/)L(z,|x')dp/,  (1) 

where  b  and  c  (in  units  of  inverse  meters)  are  the 
beam  scattering  and  attenuation  coefficients,  |3  is  the 
scattering  phase  function  normalized  such  that 
P(|x,  l)dp,  =  1,  and  p,  is  the  cosine  of  the  polar  angle 
with  respect  to  the  downward  depth  z  (in  meters). 
All  quantities  in  Eq.  (1)  are  implicitly  a  function  of 
wavelength.  The  backscattering  coefficient  bb  is  re¬ 
lated  to  b  by 


bh  =  b^  |3(|x,  l)dp.,  (2) 

and  the  absorption  coefficient  is  related  to  b  and  c  by 
a  +  b  =  c. 

The  irradiance  ratio  R{z)  and  radian ce-irradiance 
ratio  Rl(z)  are 

R(z)  =  EJzJEJz),  R\z)  =  LJz)/EJz),  (3) 

where  Lu{z )  is  the  upward  radiance  [LJz)  =  L(z, 
-1)]  and  Eu(z)  and  Ed(z )  (in  units  of  W  m-2  nm"1) 
are  the  upward  and  downward  irradiances: 

E„{z)  =  2it  |  |p-|L(z,  p.)dp,, 

EJz)  =  2tt  f  p,L(z,  p,)dp,.  (4) 

Jo 

The  upward  and  downward  diffuse  attenuation  coef¬ 
ficients  KJz)  and  KJz)  (in  inverse  meters)  are  de¬ 
fined  by 


dln[£„(z)]  dln[EJz)] 

KJz)  = - , - ,  KJz)  = - , - ,  (5) 


dz 


dz 


and  the  diffuse  attenuation  coefficient  for  L(z,  p.)  is 
denoted  by  KJz,  p.),  with 


KJz,  -1)  =  - 


d  ln[L„(z)] 
dz 


(6) 


At  large  optical  depths  in  deep  homogeneous  waters 
with  no  internal  sources,  the  values  of  R(z),  RL(z), 
and  KJz)  approach  the  asymptotic  values  R,_,  Rj ', 
and  K x,  respectively,  that  are  IOPs  of  the  water. 
Because  i?rj0,  R„L,  and  Krj0  are  IOPs,  they  can  be  de¬ 
fined  mathematically  in  terms  of  a,  b,  and  \i.  There¬ 
fore,  although  they  cannot  always  be  directly 
measured,  they  have  physical  meaning  in  any  waters 
and  are  in  general  a  function  of  depth.  At  a  partic- 
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ular  depth,  a  solution  of  the  equations  of  radiative 
transfer  yields  the  relationships3-6 

R°o  =  gj-vj/gjvj,  (7) 

R,L  =  4>(— Vi,  1  )/gJv1),  (8) 

Krj  =  c/  Vi,  (9) 

where 

gi(v)  =  |  4>(r,  pdqdp.  (10) 

Jo 

and  rq  is  the  largest  positive  eigenvalue  correspond¬ 
ing  to  the  eigenfunction  h)  of  Eq.  (1).  The  val¬ 
ues  of  iq  and  (fie,,  p)  depend  only  on  fi  and  the 
single-scattering  albedo  «0,  where  w0  =  b/c.  There¬ 
fore  Eqs.  (7)  and  (8)  are  explicit  expressions  for  R  a 
and  R./J  in  terms  of  p  and  io0  and  Eq.  (9)  is  an  explicit 
expression  for  K„  in  terms  of  |3,  b,  and  c. 

To  evaluate  v1  and  <i>(iq,  p)  from  p  and  co0,  the 
spatial  and  angular  dependencies  in  Eq.  (1)  can  be 
separated  with  the  eigenmodes  L-(z,  p)  =  <l>(± vr  p) 
exp (+cz/vj),  and  the  scattering  phase  function  can  be 
expanded  in  Legendre  polynomials,  yielding 

1  M 

P(h,  h')  =  „  2  (2n  +  1  )fnPn(\R>Prl(W),  fo  =  1,  HD 

^  n-0 

where  fn  are  the  expansion  coefficients,  Pn( p)  are  the 
Legendre  polynomials,  and  M  is  the  degree  of  scat¬ 
tering  anisotropy.  It  is  then  found  that  the  discrete 
eigenfunctions  satisfy6 

M 

4>(  ±  vj,  p)  =  °_J  2  ( 2n  +  !)/"« 

2 (vj  +  p)  „-0 

Xgn(±Vj)Pn(  p),  vj>  1.  (12) 

The  Chandrasekhar  polynomials7  gn(±v-)  obey  the 
recursion  formula 

ngjv,)  =  hn-1vJgn-Jvj)  -  (n  -  l)gn-Jvj),  (13) 

starting  withg^  =  0  and  gn  =  1,  where  hn  =  (2 n  + 
1)(1  -  oj0/„).  From  the  spherical  harmonics  (PN) 
method8  with  N  odd  and  arbitrarily  large,  the  posi¬ 
tive  eigenvalues  v7  are  approximately  the  roots  of 

&v+i(l)  =  0-  114) 


3.  Methods 

Several  steps  are  taken  to  estimate  the  values  of  a{z ) 
and  bb(z)  from  field  measurements  ofLu(z)  andErf(z). 
First,  to  reduce  noise,  Lu{z)  and  Ed[z)  are  smoothed 
with  respect  to  depth  by  making  linear  least-squares 
fits  to  In \Lu{z)\  and  ln[Erf(z)]  over  small  ranges  of 
depths.  Second,  RL(z)  and  KJz)  are  computed  from 
these  smoothed  profiles  with  Eqs.  (3)  and  (5). 

The  third  step  is  the  determination  of  RxL{z)  and 
Kx(z)  from  RL(z)  and  Kd(z).  In  optically  deep  and 
well-mixed  waters,  the  profiles  RL(z)  and  Kd{z)  as¬ 
ymptotically  approach  the  values  RXL  and  K,. 


Fig.  1.  Single-scattering  albedo  w0,  the  largest  eigenvalue  v1;  and 
the  ratio  of  scattering  to  absorption  versus  the  asymptotic 
radiance-irradiance  ratio  RjJ  for  the  Petzold  scattering  phase 
function. 


Therefore  if  it  is  known  that  the  water  is  homoge¬ 
neous  and  optically  deep,  then  Rj"{z)  and  K  Jz)  are 
independent  of  depth  and  can  be  taken  to  be  equal  to 
the  values  of  RL(z)  and  Kd(z)  deep  in  the  euphotic 
zone.  If  the  water  is  homogeneous  but  optically 
shallow  with  a  highly  absorbing  bottom,  then  R, 
and  K,  should  be  estimated  from3  RrL  maxL-R^fzjJ 
and  K,_  «<  max[iff/(z)J.  For  optically  shallow  waters 
with  a  reflective  bottom,  an  algorithm  for  the  estima¬ 
tion  of  Ry_  is  given  in  Ref.  3;  however,  an  analogous 
equation  for  R  L  could  not  be  derived  in  the  same 
manner.  In  stratified  waters,  Rj'iz)  and  K  Jz)  are 
depth  dependent  and  cannot  be  directly  measured. 
However,  Rj'iz)  and  KJz)  can  be  estimated  with 
R,J(z)  Rl{z)  and  KJz)  Kd(z).  Although  these 
seemingly  crude  approximations  are  generally  inac¬ 
curate  near  the  surface  and  are  not  adequate  to  de¬ 
termine  the  fine  vertical  structure  of  highly  stratified 
waters,  they  can  give  fairly  accurate  vertically 
smoothed  profiles  of  the  IOP’s  more  than  an  optical 
depth  from  the  surface.3 

In  the  fourth  step  the  values  of  w0(z)  and  vjz)  are 
calculated  from  Rj'(z).  In  Ref.  3,  w0  was  deter¬ 
mined  with  an  iterative  solution  method.  However, 
this  iterative  technique  is  not  necessary.  There  is  a 
one-to-one  correspondence  between  Rjg  and  co0  for  a 
given  p,  and  unique  values  of  c,  follow  from  w0. 
Therefore,  given  the  coefficients  fn  of  the  assumed 
scattering  phase  function,  a  table  can  be  generated  of 
rq  and Rj*  versus  co0  with  Eqs.  (8),  (10),  and  (12)— ( 14). 
This  table  can  be  inverted  to  obtain  w0  and  iq  versus 
R,l.  For  example,  a  plot  of  w0  and  v1  versus  R,L  for 
the  Petzold  scattering  phase  function9  is  shown  in 
Fig.  1.  Once  these  relationships  are  calculated,  w0 
and  q  can  be  quickly  determined  for  any  measured 

RooL- 

Finally,  from  vjz)  and  «0(z),  the  values  of  c(z), 
a{z),  and  b(z )  are  calculated  with  c(z)  =  vJz)KJz), 
a{z)  =  c(z)[l  —  co0(z)],  and  b{z)  =  c(z)u>0(z).  The 
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values  of  bb(z)  are  computed  from  the  substitution  of 
Eqs.  (11)  into  Eq.  (2)  to  obtain10 

M  /t0 

bb  =  (b/2)  2  (2 n  +  Dfn  |  P„(q)d|x  (15) 

=  (6/2) 


1-2  (2n  +  l)fn  |  P„(p-)dp, 

n  odd 


.  (16) 


4.  Sensitivity  Analysis 

In  practice,  the  accuracy  of  the  Lu-Ed  method  de¬ 
pends  on  the  sensitivities  of  a  and  bb  to  small  varia¬ 
tions  in  K,_  and  R  jJ.  The  normalized  sensitivity 
coefficients  of  a  and  bb  with  respect  to  K,_  are  identi¬ 
cally  unity3: 


K  ,  da  K,  dbb 
a  dK„  bh  dK„ 


(17) 


For  example,  a  1%  increase  in  K,  results  in  a  1% 
increase  in  the  calculated  values  of  a  and  bb.  Ana¬ 
lytical  expressions  for  the  sensitivity  coefficients  of  a 
and  bb  with  respect  to  R.,_,  expressed  as  functions  of 
easily  computed  forward-problem  sensitivity  coeffi¬ 
cients,  were  derived  previously3  for  the  Eu-Ed 
method.  Those  with  respect  to  Rj"  can  be  derived  in 
the  same  manner,  and  the  resulting  expressions  are 
the  same  except  with  R,_  replaced  by  R  ,L: 


R„L  da  I R„L  dv1  \  /  w0  \  / R  J'  du>0  \ 

a  dRjJ  \  iq  dR,,L)  \1  -  w0/\  w0  dRj'J  ’ 

Rj '  dbb  _  (Rj1  \  (Rj*  d«0  \ 

bb  dRjJ  \  Vi  dRjJj  +  \  w0  dRjJ]  ' 


(18) 

(19) 


The  forward-problem  sensitivity  coefficients  [( 'Rj'/ 
ViXcnq/dR./')]  and  \(RjJ/ o)0)(dw0/ dR,L)\  are  always 
positive  and  can  be  evaluated  from  the  numerical 
relationships  shown  in  Fig.  1.  It  can  be  seen  from 
Eqs.  (18)  and  (19)  that  the  effects  of  (diq/dR  ,,/J)  and 
fdoj0/clR  /J)  on  the  estimation  of  bb  are  additive, 
whereas  the  effects  of  these  on  the  estimation  of  a 
tend  to  cancel.  Therefore  a  is  much  less  sensitive  to 
R,JJ  than  is  bh.  Although  the  term  [«0/(l  —  w0  )]  in 
Eq.  (18)  is  large  for  co0  near  unity,  it  is  multiplied  by 
the  term  (dwJdR.J")  that  is  small  for  oj0  near  unity. 
Nonetheless  the  sensitivities  of  a  and  bb  to  R  ,_  and 
RjJ  are  much  larger  for  w0  near  unity  than  for  co0  < 
0.9.  The  values  of  the  sensitivity  coefficients  of  Eqs. 
(18)  and  (19)  are  shown  in  Fig.  2  for  the  Petzold 
scattering  phase  function.  Estimates  of  a  are  rela¬ 
tively  insensitive  to  R(z)  and  RL(z),  with  the  magni¬ 
tude  of  the  sensitivity  coefficients  being  less  than  0.3 
for  ojf)  <  0.9.  Although  bb  is  somewhat  more  sensi¬ 
tive  than  a  to  Rr_  and  !?/',  (R  ,Jbh)(dbh/dR,)  <  1  for 
co0  <  0.980  and  (R.,L /bb)(dbh/ dR,,L)  <  1  for  w0  < 
0.988.  Both  a  and  bb  are  slightly  less  sensitive  to  R  r_ 
than  to  R„l  for  oj0  <  0.8  and  vice  versa  for  co0  >  0.8. 
As  co0  approaches  zero,  the  normalized  sensitivity 
coefficients  of  a  and  b  approach  zero  and  unity,  re¬ 
spectively. 


Fig.  2.  Normalized  sensitivity  coefficients  ( R„/a)(da/dR „),  (R„L/ 
a)(Aa/6Rj'),  (R^/bh)(Abb/ ARJ),  and  (Rj'/bb)(dbb/dRj')  for  the 
Petzold  scattering  phase  function. 


Expressions  were  previously  determined3  for  the 
sensitivity  coefficients  of  a  and  bb  with  respect  to  a 
parameter  p  of  the  scattering  phase  function.  This 
parameter  had  been  taken  to  be  the  scattering  asym¬ 
metry  factor  but  can  be  any  scalar  descriptor  of  the 
phase  function.  The  sensitivity  coefficients  ex¬ 
pressed  in  terms  of  the  three  forward-problem  sensi¬ 
tivity  coefficients  are 

-y  =  (-  7^-V-  (20) 

a  dp  Ivi  dp  j  \1-mo/\w0  dp) 

p^dbh^p^dvi^P-dwo  +  p^dbh  ^ 

bh  dp  V!  dp  co0  dp  bb  dp  ’ 

where  du>0/dp  =  ('  1/ dR.,L)[dRrL / dp).  Estimates  of  a 
are  insensitive  to  the  scattering  phase  function,3-4  as 
the  effects  of  the  sensitivities  of  v1  and  oj0  with  re¬ 
spect  to  p  in  Eq.  (20)  tend  to  cancel.  Therefore  the 
accuracy  of  the  assumed  phase  function  is  unimpor¬ 
tant  in  the  estimation  of  a.  In  the  estimation  of  bb, 
on  the  other  hand,  the  effects  of  the  sensitivities  of  v± 
and  w0  to  p  in  Eq.  (21)  are  additive.  However,  be¬ 
cause  these  tend  to  have  the  opposite  sign  of  the  last 
term  in  Eq.  (21),  bh  is  far  less  sensitive  to  the  as¬ 
sumed  (i  than  is  b,  for  which  the  sensitivity  coefficient 
is  similar  to  Eq.  (21)  but  with  the  last  term  absent.3 
For  the  estimation  of  bb,  the  most  important  param¬ 
eter  describing  the  assumed  scattering  phase  func¬ 
tion  is  bb  =  bb/b.  If p  is  taken  to  be  bb,  then  Eq.  (21) 
becomes 


bb  dbb  _  bb  dvi  bb  du>0  . 

~  ~  ~  I  ~  I  1. 

bb  dbh  vx  dbh  co0  dbb 

Because  for  a  given  measurement  of  R,rL  an  overes¬ 
timate  of  bb  leads  to  an  underestimate  of  w0  and  vls 
and  vice  versa,  the  terms  {dvjd bb)  and  (du>0/dbb)  in 
Eq.  (22)  are  negative  and  tend  to  cancel  with  the 


20  August  1999  /  Vol.  38,  No.  24  /  APPLIED  OPTICS 


5099 


Fig.  3.  Simulated  profiles  of  the  irradiance  and  radiance- 
irradiance  ratios  and  of  the  diffuse  attenuation  coefficients  for  the 
case  of  homogeneous  optically  deep  water  with  a  =  0.30,  b  =  0.70, 
a  flat  surface,  and  an  illumination  that  is  30%  diffuse  skylight  and 
70%  direct  sunlight  at  30°  from  the  zenith. 


unity  term.  Therefore,  although  the  accurate  deter¬ 
mination  of  bb  requires  that  bh  of  the  assumed  scat¬ 
tering  phase  function  be  approximately  that  of  the 
true  scattering  phase  function,  bb  is  not  too  sensitive 
to  the  assumed  bb  to  make  its  estimation  impracti¬ 
cal.3’5 

5.  Numerical  Simulations 

Numerical  simulations  were  performed  to  evaluate 
the  Lu-Ed  method  for  determining  a(z)  and  bh(z)  and 
to  compare  this  method’s  performance  with  that  of 
the  Eu-Ed  algorithm.  Simulated  Lu{z),  Eu{z),  and 
Ed(z)  data  were  generated  for  the  San  Diego  Harbor 
Petzold  scattering  phase  function11  and  specified 
w0(z),  surface  roughness,  and  illumination  condi¬ 
tions.  Simulations  for  flat,  homogeneous  waters 
were  performed  with  the  discrete  ordinates  radiative 
transfer  code  disort,12  whereas  data  for  roughened 
surfaces  or  inhomogeneous  waters  were  generated 
with  Hydro  light. 13 

For  example,  shown  in  Fig.  3  are  ttRl(z),  R{z), 
Kd{z),  Ku(z),  and  KL(z,  — 1)  for  the  case  of  homoge¬ 
neous  optically  deep  water  with  a  =  0.30,  b  =  0.70 
( bb  =  0.014),  a  flat  surface,  and  an  illumination  that 
is  30%  diffuse  skylight  and  70%  direct  sunlight  at  30° 
from  the  zenith.  The  values  all  increased  nearly 
monotonically  with  depth.  Those  of  ttRl(z)  and  R(z) 
approached  their  respective  asymptotic  values  of 
0.011  and  0.019.  As  was  found  in  general  for  simu¬ 
lations  over  a  wide  range  of  illumination  conditions 
and  water  optical  properties,  the  profile  of  RL(z)  was 
much  more  constant  with  depth  thani?(z)  and  there¬ 
fore  closer  than  R{z)  to  its  respective  asymptotic 
value.  The  profiles  Kd(z),  Ku{z),  and  KL{z,  -1)  ap¬ 
proached  their  shared  asymptotic  value  of  0.41,  with 
Kd(z)  being  closest  to  K.a.  At  large  depths  the  values 


Fig.  4.  Estimates  of  the  absorption  and  backscattering  coeffi¬ 
cients  for  the  simulation  of  Fig.  3  obtained  with  the  Eu-Ed  and 
Lu-Ed  methods. 


of  K, o;  R  ,l,  and  R  ,_  can  be  measured  directly  with 
high  accuracy.  Therefore  if  the  waters  are  known  to 
be  homogeneous,  the  depth-independent  values  of  a 
and  bb  can  be  determined  accurately  from  deep  mea¬ 
surements  of  K,_  and  RVL  or  R.,_.  However,  if  the 
vertical  structure  of  the  waters  is  unknown,  then 
Kx(z)  and  Rj'iz)  or  R^z)  must  be  estimated  with 
KJz)  «  Kd(z)  and  Rj'iz)  «  RL(z)  or  Rx(z)  «*  R(z). 
For  the  simulation  of  Fig.  3,  the  errors  in  these  ap¬ 
proximations  at  1  m  were  —6.2%,  —2.7%,  and  —8.4%, 
respectively,  and  at  5  m  were  -0.95%,  -0.66%,  and 
- 1.3%.  The  corresponding  profiles  of  a(z)  and  bb(z) 
are  shown  in  Fig.  4.  For  the  given  conditions,  the 
retrieved  IOP  profiles  are  accurate  at  large  depths 
but  are  in  error  near  the  surface.  Estimates  of  a{z) 
from  the  Eu-Ed  method  were  slightly  better  than 
those  from  the  Lu—Ed  method.  The  errors  at  1,  3, 
and  5  m  were  -5.7%,  -2.2%,  and  -0.83%,  respec¬ 
tively,  for  the  Lu-Ed  method  and  -4.9%,  -1.9%,  and 
-0.74%  for  the  Eu-Ed  method.  Conversely,  esti¬ 
mates  of  bb  obtained  from  RL(z)  were  more  accurate 
than  those  obtained  from  R(z),  with  the  errors  at  the 
same  depths  being  —8.1%,  —3.4%,  and  —1.4%  for  the 
Lu-Ed  method  and  —12%,  —4.7%,  and  —1.9%  for  the 
Eu-Ed  method. 

Near  the  surface,  the  percent  error  in  a(z)  at  a 
given  depth  was  less  than  that  in  either  K„(z)  or 
R„(z).  This  is  because  the  errors  in  K,  and  R  ,_  were 
of  the  same  sign  near  the  surface  and,  as  demon¬ 
strated  in  Section  4,  the  sensitivity  coefficients  of  a 
with  respect  to  K,_  (unity)  and  with  respect  to  R  ,  (Fig. 
2)  are  of  opposite  signs.  Furthermore,  estimates  of  a 
near  the  surface  were  better  for  the  Eu-Ed  method 
than  for  the  Lu-Ed  method  because  the  Eu—Ed 
method  gave  less-accurate  estimates  ofi?„(z),  the  ef¬ 
fects  of  which  better  offset  those  that  are  due  to  errors 
in  K.Jz).  Conversely,  because  the  normalized  sen¬ 
sitivity  coefficients  of  bb  with  respect  to  K,  and  R  ,  are 
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both  positive,  the  errors  in  the  estimates  of  bb  were 
higher  than  in  either  K,_  and  and  estimates  of  bh 
from  the  Lu—Ed  method  were  better  than  those  from 
the  Eu-Ed  method.  At  large  depths,  on  the  other 
hand,  the  errors  in  KJz)  and  RJz)  are  due  primarily 
to  instrument  noise  and  natural  variability  and  are 
uncorrelated.  At  these  depths,  estimates  of  a  still 
tended  to  be  better  than  those  of  bb  because  a  is  less 
sensitive  to  R ^  than  is  bb,  and  the  Lu-Ed  method 
tended  to  slightly  outperform  the  Eu-Ed  method  for 
both  a  and  bb  for  co0  >  0.8  and  vice  versa  for  co0  <  0.8. 

In  general,  errors  in  a(z )  were  due  primarily  to  the 
errors  in  KJz).  Because  Kd{z)  was  used  to  estimate 
KJz)  for  both  the  Lu-Ed  and  Eu-Ed  methods,  the 
profiles  of  a(z)  obtained  from  the  two  methods  were 
similar  in  both  homogeneous  and  stratified  cases. 
Although  KJz)  =  Kr_  at  large  depths  in  homogeneous 
waters,  KJz)  should  not  be  used  to  estimate  K„  be¬ 
cause  internal  reflection  causes  the  value  of  KJz) 
near  the  surface  to  be  much  less  than  KJz)  and 
because  Eu(z )  typically  contains  more  instrument 
noise  than  EJz)  at  large  depths.  The  value  of  KJz, 
—  1)  on  the  other  hand  is  much  less  affected  by  inter¬ 
nal  reflection  and  may  be  the  preferred  estimator  of 
K,  near  the  surface  if  EJz)  is  highly  variable  because 
of  surface  waves. 

The  accuracy  of  a(z)  and  bjz)  retrieval  depended 
highly  on  the  surface  conditions,  with  the  best  results 
occurring  when  the  angular  distribution  of  the  light 
field  just  below  the  surface  was  most  near  the  asymp¬ 
totic  distribution.  For  example,  for  the  same  water 
conditions  as  for  Figs.  3  and  4  but  with  diffuse  illu¬ 
mination,  the  estimates  of  a(z)  and  bjz)  were  dra¬ 
matically  improved  over  those  shown  in  Fig.  4.  With 
the  Lu-Ed  method,  the  errors  in  a(z)  at  1  and  3  m 
were  only  —2.3%  and  0.49%  (compared  with  —5.7% 
and  -2.2%),  whereas  those  in  bb(z)  at  those  depths 
were  —3.7%  and  0.24%  (compared  with  —8.1%  and 
—3.4%).  Likewise,  somewhat  better  estimates  of 
a(z )  and  bjz)  were  obtained  when  the  sea  surface 
was  rough  or  when  the  Sun  was  low  in  the  sky  than 
when  the  surface  was  flat  or  the  Sun  was  at  local 
noon.  For  example,  for  the  same  illumination  con¬ 
dition  as  in  Fig.  4  but  with  15-m/s  winds  added, 
causing  a  roughened  surface,  the  errors  in  both  a(z) 
and  bjz)  in  the  top  5  m  were  reduced  by  approxi¬ 
mately  6%. 

Although  the  Lu-Ed  and  Eu-Ed  methods  for  the 
determination  of  a  and  bb  were  theoretically  derived 
for  homogeneous  waters,  they  can  be  applied  to  in¬ 
homogeneous  waters.  If  the  vertical  variation  of  the 
IOP’s  is  gradual,  it  may  be  possible  to  determine 
these  profiles  accurately,  whereas  if  the  vertical 
structure  is  fine,  only  depth-smoothed  profiles  can  be 
obtained.  Figure  5  shows  the  retrieved  profiles  of 
a(z )  and  bjz)  for  the  case  of  constant  a(z )  =  0.3,  b{z ) 
varying  sinusoidally  between  0.4  and  0.7,  a  flat  sur¬ 
face,  and  diffuse  illumination.  It  can  be  seen  that 
bjz)  is  retrieved  well  below  one  optical  depth.  The 
retrieved  a{z)  shows  sinusoidal  error,  but  its  maxi¬ 
mum  magnitude  is  only  approximately  3%.  The 
case  of  diffuse  illumination  is  shown  in  Fig.  5  because 


Fig.  5.  Determination  of  the  absorption  and  backscattering  coef¬ 
ficients  with  the  Lu-Ed  (solid  curve)  and  Eu-Ed  (dashed  curve) 
methods  for  a  simulation  (dotted  curve)  with  constant  a(z)  =  0.3 
and  6(z)  varying  sinusoidally  with  depth.  The  surface  was  flat 
and  the  illumination  diffuse. 


errors  that  are  due  to  surface  effects  are  minimal, 
allowing  the  effects  of  the  stratification  to  be  seen 
clearly.  For  a  sunny  sky  illumination,  the  errors  in 
the  estimated  a{z)  and  bjz)  near  the  surface  are 
similar  to  those  for  the  homogeneous  case  (e.g.,  Fig. 
4),  whereas  the  errors  at  large  depths  are  similar  to 
those  in  Fig.  5. 

6.  Field  Experiment 

Underwater  light  and  IOP  measurements  were  made 
on  6  and  8  October  1997  at  nine  stations  in  Long 
Island  Sound  between  the  mouths  of  the  Thames  and 
Connecticut  rivers.  The  locations  and  depths  of 
these  sites  are  given  in  Table  1.  The  objective  was  to 
compare  values  of  a(z)  obtained  from  the  Lu—Ed  and 
Eu-Ed  methods  with  those  from  a  reflecting  tube  in¬ 
strument.  Verification  of  bb(z)  values  from  the  light 
measurements  was  not  possible  because  there  was  no 
independent  measurement  available. 

An  upwelling  radiance  sensor  (OCR-200,  Satlantic, 
Inc.,  Halifax,  Nova  Scotia)  was  used  to  obtain  Lu{z) 
measurements,  and  a  downwelling  irradiance  sensor 
(OCI-200,  Satlantic,  Inc.)  was  flipped  vertically  be¬ 
tween  casts  to  obtain  EJz)  and  EJz)  alternatively. 
A  data-acquisition  system  (DATA-100,  Satlantic, 
Inc.)  was  used  to  determine  the  depth  and  to  merge 


Table  1.  Location  and  Depth  of  the  Cruise  Stations  in  Long  Island 
Sound 


Station  Number 

Latitude  (N) 

Longitude  (W) 

Depth  (m) 

1 

41°16.218' 

72°11.711' 

29 

2 

41°19.097' 

72°11.195' 

5 

3 

41°50.884' 

72°16.970' 

7 

4 

41°14.703' 

72°17.175' 

32 

5 

41°15.438' 

72°20.008' 

8 

6 

41°15.014' 

72°20.568' 

8 

7 

41°13.892' 

72°17.901' 

41 

8 

41°13.293' 

72°13.585' 

50 

9 

41°17.913' 

72°04.975' 

9 
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Fig.  6.  Comparison  at  site  8  of  estimates  of  the  absorption  coef¬ 
ficients  at  (from  right  to  left)  411,  443,  490,  and  555  nm  from  the 
ac-9  (solid  curve)  with  those  from  the  Lu-Ed  method  (dashed 
curve). 


the  data  streams.  At  each  site,  the  Satlantic  pack¬ 
age  was  deployed  twice,  once  to  obtain  EJz)  and 
LJz)  and  once  to  obtain  Eu{z)  and  LJz).  At  each 
wavelength,  values  of  EJz)  and  Kd(z)  were  obtained 
at  0.25-m  spacing  by  calculating  a  least-squares  lin¬ 
ear  fit  ofln[Fd(z)]  over  a  2-m  window  centered  at  each 
depth.  Likewise,  values  of Eu(z)  and  LJz)  were  de¬ 
termined  at  0.25-m  intervals  and  the  ratios  RL{z)  and 
R(z )  were  formed.  At  each  depth,  R,l,  Rr_,  and  Kx 
were  estimated  as  the  median  values  of  R{z),  RL(z), 
and  KJz),  respectively,  over  a  2-m  window  centered 
at  that  depth.  Estimates  of  a(z)  were  obtained  from 
Rj{z)  and  KJz)  with  the  methods  described  in  Sec¬ 
tion  3  and  from  R.Jz)  and  KJz)  with  the  method  of 
Ref.  3. 

Profiles  of  a(z)  were  also  obtained  with  an  ac-9 
instrument,  with  the  sensor  depth  determined  with  a 
conductivity,  temperature,  and  depth  instrument 
(Falmouth  Scientific,  Inc.,  Cataumet,  Massachu¬ 
setts).  At  each  wavelength  X,  estimates  of  ajz) 
were  made  by  first  correcting  for  temperature  ef¬ 
fects14  in  the  absorption  at  715  nm,  a7lJz),  and  then 
correcting  for  scattering  errors15  at  all  other  wave¬ 
lengths  by  subtracting  from  ajz)  the  term 
\a11Jz)bJz)/b11Jz)\.  The  contribution  of  pure  wa¬ 
ter  to  ajz)  was  added,  using  the  values  of  Pope  and 
Fry2  at  412  and  440  nm  and  those  of  Tam  and  Patel16 
for  longer  wavelengths. 

Comparisons  of  ajz)  from  the  Satlantic  measure¬ 
ments  and  from  the  ac-9  instrument  were  made  at 
the  five  shortest  Satlantic  wavelengths  (411,  443, 
490,  509,  and  555  nm).  All  stations  were  character¬ 
ized  by  calm  seas  and  clear  skies.  Four  of  the  sites 
(1,  4,  7,  and  8)  were  optically  deep  at  all  wavelengths, 
whereas  detectable  light  reached  the  bottom  at  the 
other  five  sites.  A  strong  current  made  it  impossible 
to  keep  the  Satlantic  instruments  in  a  vertical  orien¬ 
tation  at  station  5. 

Shown  in  Fig.  6  is  an  example  comparison  (from 
site  8)  of  the  estimates  of  a411(z),  a443(z),  a490(z),  and 
a555(z)  obtained  with  the  ac-9  with  those  obtained 
from  the  Lu—Ed  method.  The  ac-9  data  were  binned 


absorption  coefficient  from  Satlantic 


Fig.  7.  Absorption  coefficients  determined  with  the  Lu-Ed 
method  versus  those  obtained  with  the  ac-9.  Included  are  data  at 
2-m  intervals  from  sites  1-4  and  7-9.  Also  shown  is  the  line  aSat 

“ac-9- 


to  0.1-m  vertical  resolution.  Of  the  wavelengths 
tested,  light  at  555  nm  penetrated  the  deepest,  pro¬ 
viding  the  longest  profiles  of  LJz)  and  Ed{z).  The 
agreement  between  the  two  methods  is  good,  and  the 
gross  vertical  structure  is  provided  by  both  methods; 
however,  the  natural  light  field  approach  can  match 
neither  the  spatial  resolution  nor  the  vertical  range  of 
the  ac-9. 

Shown  in  Fig.  7  is  a  scatter  plot  of  the  values  of  the 
absorption  coefficient  obtained  with  the  Lu-Ed 
method,  aSat(z),  versus  those  obtained  with  the  ac-9, 
aac-9(z).  The  values  shown  are  for  all  stations  except 
stations  5  and  6  and  were  obtained  at  2-m  vertical 
spacing  by  taking  the  median  value  over  a  2-m  win¬ 
dow  about  each  depth.  Superimposed  in  Fig.  7  is  the 
line  of  perfect  correlation,  as,Jz)  =  aac_g(z).  A  linear 
fit  to  the  data  in  Fig.  7  yielded  a  slope  of  0.97  and  an 
offset  of  —0.007,  with  a  correlation  coefficient  of  r2  = 
0.98.  No  consistent  trend  was  determined  in  the 
ratio  aSat:aac.g  with  respect  to  depth.  Meaningful 
regressions  at  individual  wavelengths  were  not  pos¬ 
sible  because  of  the  limited  range  of  the  values  of  the 
IOP’s  within  the  data  set.  A  regression  of  the  aSat 
and  aac_9  data  for  station  5  (in  the  Connecticut  River) 
yielded  a  slope  of  1.06  and  a  correlation  coefficient  of 
r2  =  0.97;  however,  these  data  were  not  included  in 
Fig.  7  because  there  was  an  offset  of  0.29  m  ,  with 
asat  >  a  ac-9-  The  offset  was  due  presumably  to  a  tilt 
in  the  Satlantic  package  (because  of  the  current)  that 
caused  an  underestimation  of  R „  and  co0  and  thus  an 
overestimation  of  a(z).  Similarly,  the  regression  of 
the  data  at  site  6  had  a  high  value  of  r 2  but  was  offset 
by  0.08  m_1. 

The  Lu—Ed  approach  was  easier  to  implement  than 
the  Eu-Ed  method  because  the  LJz)  signal  contained 
less  instrument  noise  than  the  Eu  signal  and  the 
profiles  of  LJz)  and  EJz)  were  taken  simulta¬ 
neously,  whereas  the  EJz)  and  EJz)  profiles  were 
done  in  succession.  Nonetheless,  the  two  methods 
were  generally  in  good  agreement.  In  fact,  the  esti¬ 
mates  of  a(z)  from  the  two  methods  tracked  each 
other  closely  when  the  same  KJz)  profiles  were  used 
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Fig.  8.  Absorption  coefficients  determined  with  the  Eu-Ed 
method  versus  those  obtained  with  the  ac-9  for  sites  1,  4,  and  8. 


for  each.  Shown  in  Fig.  8  is  a  scatterplot  produced  in 
the  same  manner  as  Fig.  7  but  for  the  Eu—Ed  method 
at  sites  1,  4,  and  8.  The  regression  slope,  offset,  and 
r2  values  were  1.08,  -0.002,  and  0.98,  respectively. 

7.  Discussion 

A  method  was  proposed  for  estimating  a(z)  and  bb(z) 
from  in-water  measurements  of  Lu(z)  and  Ed(z). 
This  method  is  valuable  because  it  provides  a  large- 
volume  measurement  of  the  IOP’s  that  can  be  used  to 
obtain  optical  closure  with  small-volume  measure¬ 
ments  such  as  that  provided  by  the  ac-9  instrument. 
The  computation  of  the  IOP’s  is  fast,  requiring  no 
iteration,  and  no  information  about  the  illumination 
or  surface  conditions  is  necessary. 

This  method  is  implemented  by  first  determining 
Rxl(z)  and  K,  and  then  calculating  a(z)  and  bh(z) 
from  these  values.  A  sensitivity  analysis  shows  that 
the  estimates  of  a(z)  and  bh(z)  are  relatively  insen¬ 
sitive  to  errors  in  Rj'iz)  and  K.,  (z)  and  to  the  as¬ 
sumed  scattering  phase  function.  On  the  other 
hand,  numerical  simulations  show  that  the  accuracy 
of  estimates  of  R,  Jz)  and  K  Jz)  depends  highly  on  the 
illumination  conditions,  degree  of  vertical  stratifica¬ 
tion,  and  the  optical  depth  to  which  reliable  light 
measurements  can  be  taken. 

Optical  closure  was  attempted  with  measurements 
in  Long  Island  Sound  by  comparing  values  of  a{z) 
obtained  with  the  light  measurement  approach  with 
those  obtained  with  an  ac-9  instrument.  The  mea¬ 
surements  were  taken  on  days  with  clear  blue  skies, 
which  numerical  simulations  suggest  are  relatively 
poor  conditions  for  estimating  a(z)  near  the  surface 
from  light  measurements  (because  errors  in  the 
method  that  are  due  to  the  directional  nature  of  the 
illumination  reach  deepest  into  the  water  column  un¬ 
der  sunny  conditions).  Nonetheless,  the  agreement 
between  the  two  methods  in  the  field  was  good,  with 
estimates  of  a(z )  being  systematically  larger  than 
those  obtained  from  the  ac-9  instrument  by  only  an 
insignificant  amount.  In  the  future,  the  Lu-Ed 
method  should  be  compared  with  ac-9  measurements 


over  a  wider  range  of  IOP  values  than  those  encoun¬ 
tered  in  our  field  experiment  so  that  regressions  can 
be  made  at  individual  wavelengths. 

Application  of  the  natural  light  field  approach 
proved  to  be  difficult  at  short  wavelengths  because  of 
the  high  light  attenuation.  The  performance  of  the 
sensor  is  therefore  important  at  blue  wavelengths. 
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Stamnes  and  Curtis  Mobley,  respectively. 

References 

1.  J.  R.  V.  Zaneveld,  “A  reflecting  tube  absorption  meter,”  in 
Ocean  Optics  X,  R.  W.  Spinrad,  ed.,  Proc.  SPIE  1302, 124-136 
(1990). 

2.  R.  M.  Pope  and  E.  S.  Fry,  “Absorption  spectrum  (380-700  nm) 
of  pure  water.  II.  Integrating  cavity  measurements,”  Appl. 
Opt.  36,  8710-8723  (1997). 

3.  R.  A.  Leathers  and  N.  J.  McCormick,  “Ocean  inherent  optical 
property  estimation  from  irradiances,”  Appl.  Opt.  36,  8685- 
8698  (1997). 

4.  H.  R.  Gordon  and  G.  C.  Boynton,  “Radiance— irradiance  inver¬ 
sion  algorithm  for  estimating  the  absorption  and  backscatter- 
ing  coefficients  of  natural  waters:  homogeneous  waters,” 
Appl.  Opt.  36,  2636-2641  (1997). 

5.  H.  R.  Gordon  and  G.  C.  Boynton,  “Radiance-irradiance  inver¬ 
sion  algorithm  for  estimating  the  absorption  and  backscatter- 
ing  coefficients  of  natural  waters:  vertically  stratified  water 
bodies,”  Appl.  Opt.  37,  3886-3896  (1998). 

6.  N.  J.  McCormick,  “Analytical  transport  theory  applications  in 
optical  oceanography,”  Ann.  Nucl.  Energy  23,  381-395  (1996). 

7.  S.  Chandrasekhar,  Radiative  Transfer  (Oxford  U.  Press,  New 
York,  1950). 

8.  M.  Benassi,  R.  D.  M.  Garcia,  A.  H.  J.  Karp,  and  C.  E.  Siewert, 
“A  high-order  spherical  harmonics  solution  to  the  standard 
problem  in  radiative  transfer,”  Astrophys.  J.  280,  853-864 
(1984). 

9.  C.  D.  Mobley,  B.  Gentili,  H.  R.  Gordon,  Z.  Jin,  G.  W.  Kattawar, 
A.  Morel,  P.  Reinersman,  K.  Stamnes,  and  R.  H.  Stavn,  “Com¬ 
parison  of  numerical  models  for  computing  underwater  light 
fields,”  Appl.  Opt.  32,  7484-7504  (1993). 

10.  Z.  Tao,  N.  J.  McCormick,  and  R.  Sanchez,  “Ocean  source  and 
optical  property  estimation  using  explicit  and  implicit  algo¬ 
rithms,”  Appl.  Opt.  33,  3265-3275  (1994). 

11.  T.  J.  Petzold,  “Volume  scattering  functions  for  selected  ocean 
waters,”  SIO  Ref.  72-78  (Scripps  Institution  of  Oceanography, 
La  Jolla,  Calif.,  1972). 

12.  Z.  Jin  and  K.  Stamnes,  “Radiative  transfer  in  nonuniformly 
refracting  layered  media:  atmosphere- ocean  system,”  Appl. 
Opt.  33,  431-442  (1994). 

13.  C.  D.  Mobley,  “Hydrolight  3.1  Users’  Guide,  1996,”  final  report 
(SRI  International,  Menlo  Park,  Calif.,  1996). 

14.  W.  S.  Pegau,  D.  Gray,  and  J.  R.  V.  Zaneveld,  “Absorption  and 
attenuation  of  visible  and  near-infrared  light  in  water:  de¬ 
pendence  on  temperature  and  salinity,”  Appl.  Opt.  36,  6035- 
6046  (1997). 

15.  C.  S.  Roesler  and  J.  R.  V.  Zaneveld,  “High  resolution  vertical 
profiles  of  spectral  absorption,  attenuation,  and  scattering  co¬ 
efficients  in  highly  stratified  waters,”  in  Ocean  Optics  XII,  J.  S. 
Jaffe,  ed.,  Proc.  SPIE  2258,  309-319  (1994). 

16.  A.  C.  Tam  and  C.  K.  N.  Patel,  “Optical  absorptions  of  light  and 
heavy  water  by  laser  optoacoustic  spectroscopy,”  Appl.  Opt.  18, 
3348-3357  (1979). 


20  August  1999  /  Vol.  38,  No.  24  /  APPLIED  OPTICS 


5103 


